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Abstract 

We analyze the short and long distance contributions to inclusive B —> decay, paying 

particular attention to the dependence on the Cabibbo-Kobayashi-Maskawa parameter 
Vtd- We discuss penguin diagrams with internal u and c quarks in the framework of the 
effective held theory. We also estimate the size of possible long range contributions by 
using vector meson dominance. 
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1 Introduction 


The radiative penguin decays have recently received considerable theoretical attention. 

The short distance QCD corrections to 6 —> sy decay have been calculated completely at 
the leading order |]l| and partially at the next-to-leading order [Q ; the branching ratio is in 
agreement (within the errors) with the experimental value from CLEO 0. The penguin 
diagrams for the 6 —>■ sy decay are dominated by the virtual t-quark contribution that 
is proportional to the element of the Cabibbo-Kobayashi-Maskawa (CKM) matrix 
On the other hand, the 6 —dy decay can, if dominated by the t-quark loop, provide 
information on Vtd- This decay has not been detected yet. The expected branching ratio 
is approximately 10“^ (see e. g. |^) and the statistics required is of order 10^ — 10® B 
mesons, which is within the reach of future CLEO/B factories. The principal experimental 
challenge is to keep a good control of the main background due to the 6 —>• s y decay. 

In this paper we will analyze the inclusive B —> Xd y decay. It has been suggested (see 
e.g. p) that its amplitude may not proportional to \Vtd\ because of sizable contributions 
from penguin diagrams with internal c and u quarks. We compare the relative contribu¬ 
tions of penguin diagrams with heavy and light loops in the framework of the effective 
held theory. We are interested whether the dependence on the light masses is logarith¬ 
mic or powerlike, since we expect that terms with powerlike dependence are negligible. 

We also aim to estimate possible long range contributions to the decay rate of inclusive 
B —> X^y; we adopt the phenomenological approach of vector meson dominance (VMD). 

2 The Short Distance 

The amplitude A for the b ^ s j or b —>■ d j decay can be written as A = Vt At + A^ + Vc Ac, 
where Vt = Vth V)*, Vc = Wb Vcxi W = Wb y^xi with x = s or d respectively. By using the 
unitarity of the CKM matrix, we can rewrite the amplitude in the form 

A = Vt{At — Ac) + Vu{Au — Ac). ( 1 ) 

In the case of 6 ^ s y, we have Vt O(A^) and Vu ~ O(A^), where A ~ 0.22 is the 
Cabibbo suppression factor; thus, Vu -C Vt and it is generally considered safe to ignore the 
second term in (|^). For the b —>■ d j, however, Vt and Vu are comparable ~ 0(A®); that 
prompt us to examine more carefully the size of Au — Ac- 

Let us hrst consider the b ^ d 'y amplitude in the absence of QCD corrections. The 
gauge invariant form for Au^c,t resulting from the expansion to the second order in the 
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external momenta is 

Aq = A^ + Af q = {u,c,t}, (2) 

with 

A^ = Fp{mym^w) d{kk^ - A:V) (1 - 75)&e^ 

Ay = Ja^''fci.[mrf(l-75)+mb(l + 75)]6e^, (3) 

where is the photon polarization vector, k is the photon momentum and the coefficients 
Fp, Fm are independent of the momenta. A^ and A^ will be referred as the penguin 
and magnetic moment amplitudes. As worked out explicitly in p, the dependence of 
Ay on the ratio my my is power like, while in A^ there are also logarithmic terms. This 
can be qualitatively understood in the following way. By naive power counting in the 
euclidean space, we observe that, before the expansion in the external momenta, the 
one loop diagrams have no IR divergence because of the massive W propagator. After 
the expansion, the only source of infrared (IR) divergence in the one loop diagrams is a 
logarithmic dependence on the internal quark masses, when the masses go to zero. The 
expansion induces the following IR behaviour, expressed symbolically in terms of the 
internal momentum I 

A^-KjdH^. ( 4 ) 

Therefore we expect IR divergence (and consequently logarithms in the light masses) only 
in the penguin terms 0. When the photon is real, only the amplitudes Ay contribute to 
the decay. On the other hand. Ay and Ay are strongly suppressed due to their powerlike 
dependence on myrn^. By using the results of Inami and him ^j, we hnd 

A — A A^ — A^ 

^ Ay - A^ ^ ^ ^ 

Therefore, in the absence of QCD corrections, it is justihed to assume proportionality to 
Vt in the b ^ d j decay. 

In order to include QCD corrections, we introduce the effective held theory formalism 
and work at the lowest order in the weak interactions. The basis of operators {Qi...(58} 
for the 6 —S 7 decay is well known and we will not explicitly reproduce it here; the 
reader is referred to |]^. After the substitution of the s-quark with the d-quark, all the 

^Note that this is not necessarily true in diagrams with more than one loop, since the presence of 
more internal loop momenta and of IR subdivergences alter the naive power counting of (^). 
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operators in this basis become suitable for the 6 —c? 7 decay. In addition, we have two 
new current-current operators generated by integrating out the hh-boson 

Ql = da7t,i^ - lb)Ul3 h/37^(l - 75)6a 

- 'y5)ua h/ 37 /,(l - 75 ) 6 / 3 . ( 6 ) 

Since the effective theory and the full theory share, by definition, the same low energy 
behaviour, all possible logarithmic singularity in the light masses (for rric, mu 0 ) are 
cancelled at the matching. Therefore, at the matching scale p = M^, the coefficients 
of the effective hamiltonian may have a powerlike and logarithmic dependence on the 
heavy masses, but only a powerlike dependence on the light masses. The only non-zero 
coefficients are 

C2{mw) = 

Cj{mw) = 


C^{mw) — 


C^mw) = 

where F 2 and F 2 are Inami-Lim coefficients 0. The QCD rescaling does not change the 
powerlike dependence of the amplitude on the light masses. The anomalous dimension 
matrix for 6 ^ d 7 differs from the anomalous dimension matrix for 6 —»• s 7 only for 
the additional entries due to and Q^- Thanks to their similarity with Qi and Q 2 , Qi 
and Q 2 mix with all the other operators in the same way. Due to the flavor difference, 
they do not mix with Qi and Q 2 - In 6 —>■ s 7 and b ^ d '-f, the amplitude is proportional 
to the so-called effective coefficient Cy® (see e.g. | 0 ]), that includes finite contributions 
from matrix elements. At the leading order and in naive dimensional regularization, the 
only non-zero one-loop matrix elements are the matrix elements of Qa, <^ 4 , Q 5 , Qe with a 
massive internal 6 -quark. These finite terms have been calculated in the b ^ s j case 

^ F2 is given by Eq. (B. 3 ) in |^. F2 comes from the diagrams (a), ( 6 ), (c) and (d) in Q when the Z 
is replaced by a gluon; therefore F2 = QF2 + .... 
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and are obviously left invariant by the inclusion of and Q 2 in the basis. There are no 
other possible sources of logarithms in the light masses at the leading order; as a result, 


IS 


Cf(p) = C,{mw) 

{C2{mw) + C 2 (mw)) X! 


+ 


( 8 ) 


where rj = agirny /)/The numbers a* and are given by the eigenvectors and 
eigenvalues of the anomalous dimension matrix |^. By equation (^) and using the uni- 
tarity of the CKM matrix, we can see that the expression is still proportional to Vt, as 
pointed out also in [§. 

It must be stressed that this result is peculiar to 6 —>• d 7 (and to 6 —>• s 7 ). For 
instance, in decays like 6 —> {d,s) qq there are non-zero penguin diagrams with a virtual 
photon and massless internal quarks. They have a logarithmic dependence on the light 
masses and introduce non-negligible contributions proportional to Vu in the amplitude . 

At the next-to-leading order the situation is different. Now the matching is performed 
at as order with two loops diagrams. The expansion in the external momenta of the 
two loops diagrams gives rise to penguin and magnetic moment amplitudes like in (§) 
and d^). Since the naive power counting in (H) does not necessarily work at two loops, 
both amplitudes may have a logarithmic dependence on the internal masses. We expect 
the singularity in the light masses to be reproduced in the effective theory by the matrix 
elements of the four quark operators. This is similar to what happens in 6 —> s e“; the 
logarithms cancelled at the matching in the coefficient of the operator 57 ^( 1 — 75)6 e'y^e are 
recovered at lower energy in the matrix elements of Hef f H. At the same way, in 6 —> d 7 , 
nothing prevents the two loop matrix elements of Qi, Q 2 , Qi, Q 2 (proportional to Vc and 
Vu) to have a non-negligible logarithmic dependence on the light masses. Unfortunately, 
the next-to-leading calculation has not been completed even for the b —>■ s j decay, 
making a reliable estimate difficult; roughly speaking, we expect non-negligible terms 
proportional to Vu in (|l]) to be suppressed by an extra power of as- Since as{mb) ~ 0 . 21 , 
however, this contribution may well be sizable. 

In the inclusive B 7 decay one has to take into account also the gluon bremsstrahlung 

corrections via the b ^ d + g + j decay. These corrections have been partially calculated 
. The matrix elements of the four fermion operators in b^d + g + 'j introduce 


m 


additional terms that are not proportional to Vtd with a logarithmic dependence on the 
light masses. 
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3 The Long Distance 


The short distance approach leaves out possible contributions due to intermediate hadronic 
states. These contributions may be estimated by using the vector meson dominance 
(VMD) hypothesis, as suggested by |]^. According to VMD, the b ^ dj decay is de¬ 


scribed by the b ^ dV decay, where Id is a vector meson, (-^ and its excited states, p and 
a;), followed by the conversion Id —> 7 . Since the applicability itself of VMD is not well 
established in the case of radiative penguin decays, we will limit ourselves to give an order 
of magnitude of this contribution and we will make very simple hypothesis. By Lorentz 
invariance, the most general interaction proportional to the vector meson polarization 
is given by 


A oc ej'd{ 7 ^[ai ( 1 -75)-h&i (1-^75)] 

cr^''[a2 (1 - 75) 62 (1 + 75 )JPiy + (1 - bs) + &3 (1 + 75 )]qu 

+ [04 (1 — 75) -|- 64 (1 -|- 75)]^^ + [05 (1 ~ 75) + ^5 (1 + 75)]*?^}^) ( 9 ) 

where Pfi = pfi+ p^^ is the sum of d and b momenta and Qfi = pf^— p^^ is the vector meson 
moment. Terms proportional to sq^{l ± 75 ) b are zero since the vector meson is on-shell 
and = 0. Since the quarks are on shell, we may use the Gordon decomposition to 
write 


A oc e^d{a^'^[a2 (1 — 75) -|- 62 (1 + 75 )]Pu + (1 ~ 75) + ^3 (1 + 75)]?;^ 

-|- [04 (1 — 75) -|- 64 (1 -|- 75)] p ^}&- 


( 10 ) 


The gauge invariance requires that only transverse terms couple to the photon; therefore, 
we obtain the following amplitude 


Aj'ib 


Gf 

V2‘ 


dV) = ^Vqfvimv 


my 
rub L 


as da^^^l - 75 ) -h 63 da^''{l + 75 ) 




( 11 ) 


where fv is the Id decay constant (see Table 1) and ey is the charge factor of the con¬ 
stituent states according to the quark model (cp = l/\/2, e^j = 1/(3 a / 2 ), ep = 2/3). 
The two terms in (|TI|) do not interfere (m^ ~ 0 ) and give the same contribution to the 
transverse rate: |A'r(& —> 57 )/ oc where ^y = \J\a^\^ + I 63 /. 

Using the measured transverse decay rate for B ^ ip X, we can fit Experimen¬ 
tally [ 0 , [g 


Br{B^ip + X) = 0.81 ±0.08% 
Tl/T = 0.78 ±0.17 
Tl/T = 0.59 ±0.15 
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(1.4 GeV <pp< 2.0 GeV), 
{pp < 2.0 GeV), 


( 12 ) 

(13) 

(14) 







where is the longitudinal decay rate. The branching ratio (|^ is for direct production 


of 'll), while the polarization data (M) and (1^) include decay products of heavier charmo- 
niums. It is known that the momentum range 1.4 GeV < < 2.0 GeV is dominated 

by the exclusive modes ip K{K*), while the low momentum region < 1.1 GeV is mostly 
from heavier charmoniums. We expect (|T^) to be an overestimate since ipK is purely 
longitudinal by kinematics, and ( 0 ) to be an underestimate since it is diluted by decay 
products. Here, we simply take the average and use F^/F = 0.7 for direct inclusive ip 
production. Using Eqs (pd]), and Fy/F = 1 — F^/F = 0.3 we obtain ~ 0.19. 

The measured branching ratio for the inclusive B decay into ip' is Br(B ip' + X) = 
(0.34 ± 0.04 ± 0.03) % 1]^. We assume that the measured branching ratio is equal to 
the direct branching ratio, since there is no known cascade process for ip' production. By 
using again Ft/F = 0.3, we estimate ~ 0.18. Experimental data are not available yet 
for the other ip resonances, p and ui. Encouraged by and above, we assume the 
same as, 63 and take iv = ip = y^jasp + \h^\^ ~ 0.2 for all vector mesons V = ip, ....p, uj. 
Then the VMD amplitude for the charmonium states is 


Aip{h 


dl)= 




'fm 


as - 75) Qub + 63 + 75) q^b 


M' 


(15) 

Note that now the decay constants are extended to q‘^ 0. That gives rise to a suppression 

factor /|.(0) = where we take k ~ 0.12 from Ref. [^, no suppression 

factor is taken for p and oj, whose masses are smaller. By substituting in ([T^) Vc —»• v^, 
ep ep(^)l/\/2, /|(m%) ^ /p(^)(”^p), k ^ 1 we obtain A^^‘^\b ^ d'j). 

The contributions to the amplitude coming from ip and its resonances are approx¬ 
imately equal to the contributions coming from p and u) intermediate states times the 
GKM factors 

\^T + ^tL, 1 n-11 n . 


|h: 


~ 1.06 


~ 0.4. 


The sign between the ip mode and the uj + p mode may not be reliable, since, at least in 
the exclusive case ip K*, we expect sizable hnal state interactions. However, if we ignore 
possible hnal state phases, the contribution proportional to Vu coming from charmonium 
states (by unitarity Vc = —Vt — v^) approximately cancel the contribution coming from p 
and u) [|ll| . 

The short distance amplitude at the leading order in QGD corrections is 


G 


1 


A^^{b ^dj)=e + 75) q^bel. 


V2 


47r^ 


( 17 ) 
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where {mi,) ~ 0.3 from Eqs (|^ and (|). 

The decay rate from ([T^) alone is abont 10“^ of the short distance decay rate, in 
agreement with []^. The change to decay rate from short distance due to long distance 
contributions proportional to is less than 6% with respect to the short distance decay 
rate, while the change due to long distance contributions proportional to Vu is less than 
12% X |n„|/|nt|. The change in the decay rate due to A(^ + At‘^ is less than 6% x |n.u|/|nt|. 
Therefore, the ratio between the short and the long distance pieces is small with respect 
to the theoretical errors in the short distance decay rate itself which has been estimated 
to be about 25% [fitl. 
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V 

Tv^e+e- (GeV) lli 

fv (GeV) 

(0.770) 

6.74 X 10"^ 

0.216 

0^(0.782) 

0.60 X 10-^^ 

0.194 

i/>(3.097) 

5.27 X 10-^^ 

0.406 

ij{3.685) 

2.44 X 10-® 

0.301 

V’(3.770) 

2.83 X 10"^ 

0.104 

V’(4.040) 

7.28 X 10"^ 

0.172 

V’(4.160) 

0.77 X 10-« 

0.180 

?/>(4.415) 

0.47 X 10"^^ 

0.145 


Table 1: Values of the measured decay rate of the vector meson V into e~ pair and of the 
decay constant fy. The V decay constant is dehned by < Q\J^\V >= ey my fv{my)e^ 
and is calculated from the experimental data for r(V —> e"*" e“), according to the formula 
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